Quantum correlations from Brownian diffusion of chaotic level-spacings 
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Quantum chaos is linked to Brownian diffusion of the underlying quantum energy level-spacing 
sequences. The level-spacings viewed as functions of their order execute random walks which imply 
uncorrelated random increments of the level-spacings while the integrability to chaos transition 
becomes a change from Poisson to Gauss statistics for the level-spacing increments. This universal 
nature of quantum chaotic spectral correlations is numerically demonstrated for eigenvalues from 
random tight binding lattices and for zeros of the Riemann zeta function. 
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Over the last couple of decades a general consensus 
has emerged concerning remnants of classical chaos in 
quantum physics. Since the smooth and reversible quan- 
tum evolution allows no chaos what is nowadays called 
quantum chaos had to be traced elsewhere, usually in the 
stationary eigenvalues of a quantum system 1]. A key 
quantity which was introduced to distinguish between 
quantum chaos and quantum integrability is the distri- 
bution function of level-spacings between nearest eigen- 
values 0, 0] ■ The quantum chaotic behavior is then de- 
fined as the resemblance of the level-spacing statistics to 
those obtained from Gaussian random matrices having 
highly correlated levels, which are described by random 
matrix theory (RMT) [111110. The major character- 
istic of quantum chaos is the so-called "level-repulsion" 
between energy levels which is due to the dramatic low- 
ering of symmetries in chaotic systems, while the corre- 
sponding level-spacing distribution obeys the Wigner sur- 
mise. On the contrary, quantum systems with integrable 
classical analogues are described by Poisson level-spacing 
distribution function which implies uncorrelated random 
level-spacings with possible "degeneracies" related to the 
presence of symmetries. Quantum chaos has been ex- 
perimentally observed for the irregular energy spectra in 
atom-optics billiards HQ- 

The present era of quantum computation, where the 
exploitation of quantum information resources came to 
the forefront, requires a better understanding of spec- 
tral noise arising from quantum chaotic systems. This 
noise could be related to quantum correlations, such as 
entanglement which is measured via von Neumann en- 
tropy and usually violates Bell inequalities jl, |1|. For 
such a purpose the energy level-spacing statistical distri- 
bution function is insufficient to describe fully a quan- 
tum system. In order to reveal more about the nature of 
quantum spectral noise one must examine level-to-level 
correlations. This was done in where they conjec- 
ture that quantum chaotic energy spectra are character- 



ized by l/ f noise, which is also verified by numerical 
data [lOj . In this paper we do not look only correlations 
between the enrgy eigenvalues but also view the nearest 
level-spacings as functions of their order, so that they can 
imitate discrete time series with their order correspond- 
ing to time. We provide evidence that in quantum chaos 
the classically chaotic temporal evolution is replaced by 
diffusion of the level-spacing dynamics. In other words, 
we can interpret the evolution of level-spacings (not the 
energies themselves) as that of fictitious particles execut- 
ing random walks in one-dimensional space. Moreover, 
we arrive at a more precise definition of quantum chaos 
which includes quantum correlations in addition to sta- 
tistical behavior. 
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The energy level series of a quantum system make up 
an irregular continuous but non-differentiable self-affine 
fractal curve described by a fractal dimension D. For 
quantum chaotic levels the noisy curve is space filling 
with D ~ 2 while for energy levels from quantum in- 
tegrable systems the absence of correlations implies less 
irregular curves with fractal dimension D = 1.5. The 
corresponding power-spectra P{f) oc /~" are character- 
ized by a power-law exponent a which obeys the dimen- 
sional relation D = {5 — a)/2 |11|. Therefore, quan- 
tum chaos implies l/f spectral noise with a = 1 while 
quantum i nteg rability gives 1//^ noise with a = 2, re- 
spectively [lOj. In order to provide a better understand- 
ing for the nature of correlations in quantum chaotic en- 
ergy spectra we ask the question: "What is the nature 
of correlations when dealing with the level-spacings, in- 
stead?" . Our main finding is that the correlations in the 
level-spacing sequences of quantum chaotic spectra are 
Brownian. Nevertheless, these are quantum correlations 
which occur in addition to the unpredictability associ- 
ated to the classical chaotic behavior. Therefore, a link 
can be established between quantum chaos and quan- 
tum correlations such as entanglement 0. These facts 
are illustrated for two basic examples of quantum chaotic 
behavior, the weakly disordered lattice and the Riemann 
zeros. Furthermore, quantum chaos is related to a Gaus- 
sian distribution for the level- spacing increments while 
quantum integrability is tight up to a Poisson distribu- 
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FIG. 1: (Color online) Energy-level sequences e„ with sub- 
tracted mean spacings for chaotic, integrable and quasi- 
integrable systems obtained from disordered and quasi- 
periodic lattices. The decreasing irregularity of the curves 
can be seen with the space filling chaotic curve approaching 
D = 2, the integrable curve approaching D — 1.5 and D close 
to 1 for the chosen quasi-integrable model. The dimensions 
shown in the figure differ from the expected values since they 
were obtained by box-counting algorithms which are usually 
not very accurate for such self-affine curves. Quantum corre- 
lations exist in (a) for the maximally noisy quantum chaotic 
curves. The quantum correlations die out in (b),(c). 



tion, respectively. This allows a simple description of the 
integrability to chaos transition as a change of the level- 
spacing increments from the discrete rare event Poisson 
distribution to the continuous Gauss distribution. 

In Fig. 1 energy level signals from electrons in lattices 
as a function of the level-index n are displayed. The nth 
level is expressed as a sum of level-spacings, beginning 
from the left end of the spectrum up to the rtth level, via 



From each level-spacing Si we have subtracted a local 
mean level-spacing (S) by following the usual unfolding 
procedure for the density of states to become constant on 
average everywhere in the spectrum. The chaotic, quasi- 
integrable and integrable energy-level series can be seen 
to display very different curve roughness. Very rough 
curves with D approaching 2 are obtained for chaotic 
levels, the roughness is smaller with D close to 1.5 for 
integrable ones while it becomes much lower, closer to 
_D = 1, for the studied quasi-integrable system where 
the corresponding curve approaches a line. In Fig. 2 
the corresponding Fourier analysis of these energy-level 
signals is shown via the power-spectrum 



P{f) 



1 

N 



N 

E 

n=l 



£„ exp 



.27r/n 



ntegrable disordered 




10' 



0(0 



10 



f 



100 



1000 



FIG. 2: (Color online) The averaged power-spectra of the 
energy levels for weakly disordered three-dimensional lattice 
systems, for a quasi-integrable critical model and for the spec- 
trum of zeros of the Riemann zeta function. 



averaged over many configurations. We show the energy 
levels from weakly disordered solids and Riemann zeroes 
which give 1// noise. Our findings for the presence of 
1// power-spectra for quantum chaotic levels agree with 
the conjecture of Ref. We emphasize that quantum 
correlations exist only for the maximally noisy spectrum 
of quantum chaotic levels of Fig. 1(a). The integrable 
system of Fig. 1(c) does not show any correlations de- 
scribed by a less rough curve with Brownian l//^noise. 
The quasi-integrable system of Fig. 1(b) does not display 
any quantum correlations either. Below we shall demon- 
strate that Brownian noise describes quantum correla- 
tions for chaotic systems not for the energy levels but for 
the level-spacing sequences, instead. 

The 1// noise found to describe the quantum chaotic 
levels is a common law for many classical complex phe- 
nomena where very rare intense effects require the coop- 
eration of various favorable cases. It can be also viewed 
as a scaling law since a change of the intensity combined 
with contraction in time leaves its form unchanged. The 
1/f law is very old and appears in many cases rang- 
ing from fluid flows, stock market price indices, earth- 
quakes, etc. However, despite of its generality its 
origin remains a mystery without any widely accepted 
general mechanism for its interpretation. The most com- 
mon explanation of 1/f noise relies on the extension of 
Lorentzian correlation function j/{uj'^ + 7^) from usual 
random telegraph noise of a single fluctuator with expo- 
nential relaxation 7, to several fluctuators with a broad 
distribution of relaxation rates. Although hard to jus- 
tify physically if ^05 (7) is uniformly distributed the aver- 
age over Lorentzians with the chosen broad distribution 
of 7's gives 1/f noise Another mechanism of 1/f 
noise seems more appropriate for quantum chaotic lev- 
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FIG. 3: (Color online) Level-spacing sequences Sn with unit 
mean-level spacing (5) — 1 for chaotic, quasi-integrable and 
integrable systems. On the left the corresponding normalized 
level-spacing distribution function P{S) can be seen. 



els. It was proposed (see for example Ref. 0|) that a 
possible origin of l/f noise can be Brownian correlated 
interevents between narrow pulses. Therefore, in order to 
obtain l/f noise for the chaotic levels e„ the correspond- 
ing level-spacings Sn which appear in Eq. (1) should 
follow a Brownian random walk and the level- spacing 
increments (spacings of spacings) S„ — Sn-i to be ran- 
dom uncorrelated having white type of noise. This can be 
understood rather easily since Brownian 1//^ noise of the 
levels arises from uncorrelated /° white noise of the level- 
spacings, similarly l/f noise observed for chaotic levels 
can arise from uncorrelated white noise level-spacing in- 
crements. 

Our main results for the level -spacing series S^, — 
En — £n-i are displayed in Fig. 3 together with the cor- 
responding level-spacing distributions P{S) on the left 
of the figure. It is well-known that from the chaotic se- 
ries the obtained ^(5*) obeys the Wigner surmise, the 
less irregular integrable series is described by a Pois- 
son distribution, while the quasi-integrable ones follow 
an intermediate semi-Poisson distribution |1!tI |. The dif- 
ference between the amount of noise displayed in the var- 
ious cases is immediately obvious. In order to show that 
quantum correlations are included in the chaotic level- 
spacings which follow Brownian motion we have looked 
at the statistical distribution of the level-spacing differ- 
ences Sn — Sn-1- For chaotic levels it should be Gaussian 
and for integrable levels of Poisson type which correspond 
to correlated and uncorrelated levels, respectively. The 
two behaviors are clearly distinguished in Fig. 4. We 
have verified this picture by computing all the levels, the 
level-spacings and the increments of level-spacings in lat- 
tice models and the Riemann zeros. In Fig. 4 we show 
the statistics of the level-spacing increments Sn — Sn~i- 
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FIG. 4: (Color online) The distribution of the level-spacing in- 
crements. The transition from integrability to quantum chaos 
becomes a transition from Poisson to Gauss. For weak disor- 
der the obtained variance of the Gaussian is ~ 0.71 and for 
the Riemann zeros is ~ 0.48. The quasi-integrable system 
shown in the inset is described by the simple Poisson e~^^ . 



For weakly disordered lattices where (3=1 and 30000 
Riemann zeroes where (3 = 2 the data can be de- 
scribed by simple Gauss distributions for both cases. In 
the integrable case the statistics is Poisson while a sim- 
ilar result of different Poisson type is obtained for the 
quasi-integrable statistics in the inset. 

We summarize our numerical computations of eigenval- 
ues responsible for the level-spacings. First, we consider 
the weakly disordered Anderson model studied on 
sites of a 3D lattice. The diagonal matrix elements of the 
corresponding are random taken from a uniform prob- 
ability distribution of zero mean and variance W'^/12 
and the off-diagonal matrix elements connecting near- 
est neighbors are unity. The short ranged and sparse 
real symmetric matrices are diagonalized to obtain long 
enough signals of raw eigenvalues En ■ For such purpose 
the Lanczos algorithm was employed and the scaling be- 
havior is achieved by considering linear lattice size up to 
L = 25 for many random configurations. In the unfold- 
ing procedure one replaces each level-spacing En+i — En 
by En+i - En = M{En+i) - N{En), where e„ = M{En) 
are the unfolded levels, with N{E) the averaged inte- 
grated density of states. In this procedure the mean level- 
spacing becomes {S) = 1 throughout the spectrum. We 
have taken linear size L = 25 for the chaotic levels disor- 
der W = IQ for 57 random configurations total of 147000 
eigenvalues, for the critical levels disorder Wc = 16.5 
for 80 random configurations total of 140000 eigenvalues 
and the integrable case disorder W = total of 10763 
random configurations total of 102800 eigenvalues. The 
averaged power spectra which are shown in Fig. 2 for an 
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energy window [0, 2] which include approximately 2000 
eigenvalues. We have estimated the power spectrum laws 
a = 2.0 ± 0.07 for the chaotic case, a = 1.52 ± 0.04 for 
the integrable case while a law close to a = 1.62 ± 0.08 
was found for the critical case (not shown). 

We have also diagonalized a simpler critical model re- 
placing the critical case of 3D disordered systems to con- 
sider the behavior at the borderline from integrability 
to chaos. This is the one-dimensional quasi-integrable 
Harper model where the site potential is Vn = 
2cos(27rCTn) with irrational a = {Vb — l)/2. This model 
displays critical behavior with multifractal spectra and 
wave functions. We find that the averaged power spectra 
gave power-law estimates close to a = 3.3 to 3.8 as the 
system size increased. The level-spacings of adjacent Rie- 
mann complex zeroes [l^ also exhibit fluctuations simi- 
lar to the ones of the energy spectra of electrons in dis- 
ordered lattices. The Riemann hypothesis states that 
all non-trivial zeros of the Riemann C(s) = func- 
tion have real part equal to 1/2 with (n = 1/2 -I- e„. 
Although numerical computations agree with Riemann 
hypothesis a full theoretical proof is not known. How- 
ever, the RMT is believed to be valid in this case so that 
the spectrum of Riemann zeros represents a hypothetical 
quantum chaotic system in the unitary universality class 
described by GUE, which corresponds to broken time- 
reversal invariance via a magnetic field. The P{S) statis- 
tics for small S in this case is oc with /3 = 2, compared 
to P = 1 for the orthogonal case of disordered systems. 
In RMT language we have a crossover from GOE to GUE 
statistics. In Fig. 2 we demonstrate 1// noise also for 
sequences of the Riemann zeroes. We summarize that 
our findings for fractal fluctuations and 1// noise for the 
chaotic spectra are linked to quantum fluctuations. 

It is reasonable to expect that fingerprints of chaos 
present in the classical time evolution can be also found 
in the stationary energy levels of the corresponding quan- 
tum system. In addition to classical chaos the chaotic 
spectra should also contain some form of quantum corre- 
lations. Therefore, it is worth searching for the two ingre- 
dients of quantum chaos: "correlations" for the quantum 
part and "randomness" for the chaotic part. We find 
that both appear in the self-affine fractal noise of the 
spectral curves which is exemplified via the random walk 
behavior found for the level-spacings. Although quan- 
tum evolution is very different from that of classical point 
particles in phase space Brownian diffusion also occurs in 
the quantum case. In Hilbert space level-spacing evolu- 
tion is important, which seems to play the same role for 
quantum chaos as the time evolution in phase space for 
classical chaos. Moreover, quantum chaos implies ran- 
dom uncorrelated level-spacing increments which obey a 
Gaussian distribution. In other words, quantum mem- 
ory effects in the evolution, such as dependence on the 
initial wave packet superpositions, eventually die out for 



the level-spacing increments where quantum chaos meets 
classical chaos. The quantum chaotic level-spacing incre- 
ments obey Gauss distribution distinguished from inte- 
grable ones which obey Poisson distribution. 

We have resolved numerically important issues con- 
cerning the nature of spectral correlations and noise in 
quantum chaos. Common signatures of quantum chaos 
such as "level-repulsion" and "spectral rigidity" arise 
from quantum "correlations" which appear in chaotic 
level-spacings. If combined with "randomness" , which 
might have a classical chaotic origin, make up what we 
call quantum chaos. On the other hand, integrability 
usually concerns uncorrelated random levels where the 
level-spacings obeys the Poisson distribution. We should 
add that integrability occurs either for pure ballistic sys- 
tems where quantum correlations exist without random- 
ness or strongly disordered systems where due to local- 
ization in random positions randomness exists without 
correlations since localized wave functions do not commu- 
nicate. Quantum chaos requires both ingredients, "quan- 
tum correlations" and "randomness" , simultaneously, for 
example both of them are present in RMT spectra. 

In summary, starting from the identification of pro- 
posed 1// spectral fluctuations in quantum chaos numer- 
ical evidence was given in favor of random walk behavior 
for the level-spacings. We also argued that the defini- 
tion of quantum chaos requires both "quantum correla- 
tions" and "randomness" . Our numerical study suggests 
a connection between chaotic motion of classical parti- 
cles in phase space and motion of the level-spacings be- 
tween eigenvalues in Hilbert space. In the quantum case 
Brownian correlations are shown to contribute to quan- 
tum effects. Although the level-spacing distribution is 
widely accepted as the best diagnostic tool of quantum 
chaos there is no direct theory attached to it, as it exists 
for example for other level correlation functions via RMT 
(see however the recent study of jl3|)- This might explain 
why the simple random walk behavior of level-spacings 
in quantum chaotic systems, where quantum correlations 
and unpredictable behavior are combined was not paid 
much attention before. 
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